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ABSIERAGT 


A aigital computer code was developed to simulate the 
transient hehaviour of a steam generator typically used in 
pressurized water reactor nuclear power plants, 

The three fundamental equations of conservation of 
mass , energy and momentum were converted to a set of linear 
algebraic equatiois using finite difference method, from 
which an explicit numerical scheme was formulted. 

On the primary side single phase flow and forced 
convection heat transfer were assumed and on the secondary 
side the possibility of the existence of different flow 
regimes was assumed and the heat transfer in these regimes 
was calculated using appropriate correlations. 



CHjiPTER 1 


INTRODUCTION 


1.1 Introducticsn.; 

With the present interest in the safety studies of nuclear 
reactor behaviour, it becomes important to study the behaviour 
of various reactor components. Steam generator forms an integral 
part o£ the nuclear power plants where heat is exchanged between 
the primary coolant, which flows through the core of the nuclear 
reactor, and the secondary fluid, which provides power for the 
turbine. The steam generator presents complexity in modeling 
because, even under normal conditions both subcooled and boiling 
heat transfer occur on the secondary side. 

Modeling of a large system via state variable methods is 
quite difficult because of "the fact that the overall system 
has to be broken down into separate models for each of the 
sub-systems and then these sub-sysbem models have to be coupled 
to construct the overall system model. 

Dynamic modeling plays an important role particularly in 
nuclear power plants because of stringent safety requirements 
to ensure failsafe operation of the plant, Thou^ a great deal 
of care io taken to ensure public safety in nuclear power plants 
there is always a finite probability that a system may fail and 
lead to an undesirable situation. Transient analysis allows us 
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to stady tlie degree of danger one may encounter due to the 
failure of an equipment under a disturbed condition and it 
gives guidelines for proper instrumentation to prevent such 
accidents. In case of a steam generator a large number of 
possibilities which may perturb its steady operation can be 
enumerated. The steam generator has primary and secondary 
circuit. In the primary circuit there may be reduced flow 
rate of liquid due to loss of coolant accident (lOCA) wh dh 
generally occurs due to the rupture in the pipe, discharge 
from the high-prossure system to a low-pressure system. There 
may also be loss of flow accidents (LOPA) due to the improper 
opening of valves and pump failure . This disturbance in coo- 
lant flow in turn affects the pressure and the temperature of 
of the liquid in the primary heat transport system. Similarly 
the change in the feed water flow in the secondary circuit 
will alter the steam output and this may upset the operation of 
the turbines which will affect the power output of the whole 
power station. So transient analysis of ‘the steam generator 
is essential to ensure that it can accommodate a certain amount 
of perturbation and does not fail to complicate the situation 
if there is an accident. 

1,2 Review of Dynamic Models for Steam G-enerators; 

In the literature there are no simplified models available 
for simulating the steam generators excepting the complicated 
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models of the reactor thermal hydraulic computer codes such 
as RBLAP -4 [l] and PLiiSH [ 2 ] etc, 

Recently a mathematical model was developed by irwood and 
Kerim [3] to simulate the U-tube steam generator. This paper 
uses linearized equatims to predict the behaviour around a 
steady state operating point. The detailed heat transfer 
regimes are not considered explicitly. 

Hoeld [4] developed a nonlinear transient model for the 
calculation of the dynamic behaviour of a vertical natural- 
circulation U-tube steam generator (UTSG-) ♦ But modern nuclear 
power plants use IBUTSCr which is an improved version of UTSG 
and has higher heat transfer effectiveness. 

I.M.D. Silva [5] developed a model for simulating the 
transients of an IBUTSGr. This model takes various flow regimes 
BUd different heat transfer regimes into accoxmt and follows 
an implicit numerical technique for calculating the dynamic 
behaviour of the steam generator. 

1.3 Present Work; 

The present work consists of developing a model to simulate 
the transient behaviour of an integral economizer U-tube steam 
generator (IBUTSG) utilized in a 600 MW^ PWR. To give an 
accurate representation of the steam generator, due considerations 
have been given to the aspects of two-phase flow and heat transfer. 
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As fundamental equations the well-kno-wn conservation 
equations for mass , energy and momentum have teen chosen. 
Additional approximation equations for the thermo d 3 mamic 
properties of water and steam, for a two-phase correlation, 
for one and two-phase friction coefficients and for the heat 
conduction coefficients have "been established. 

Por steady state (starting) values the model and computer 
code developed by Armando Oosta Pinto [6] has been used. 


CHAPTBR 2 


THE STEM GSmATQR 

2#1 A brief descripticai of the steam generator [7]: 

A schematic of the steam generator rs given in figure 2*1* 

Various components that constitute the steam generator arei 

( a) Primary Side : 

The hot primary fluid coining from the nuclear reactor 
enters the steam generator at the primary inlet plenum. 

The fluid then enters the upward portion (hot leg) of the 
TT-tubes. The primary flow continues into the downward 
side (cold leg) and into the economizer section. 

(b) Tube Metal Walls: 

The tube metal walls separate the primary and the 
secondary fluids. Heat is transferred from the primary 
to the secondary fluid through the tube walls at a rate 
determined by the overall heat transfer coefficient and 
the temperature difference. 

( c) Secondary Side; 

In this side there is an economizer just above the 
bottom plenum. This economizer is nothing but a closed 
container which cavers a portion of the cold legs and 
separates them from the rest of the secondary side. The 




Fia.2-1 A schematic of an Integral Economi 
U-Tube Steam Generator (IE UTSG) 
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feedwater enters the evaporator throu.^ this economizer, 

fhe IBUTSG design is a modification of the recirculation 
type U~tube steam generator (UTSG) used in many operating PWa's. 
The design change has been made to increase the heat transfer 
effectiveness of the system. The standard UTSGf- design is similar 
to that of the ISUTSG-, but all of the feedwater is introduced 
into the downcomer and there is no eccnomizer section. Because 
of the hi^ recirculation ratio, the water entering the heated 
section from the doxmcomer is near the saturation temperatxxre. 

Aa a result, there is a small temperature difference between 
the primary water of the cold legs and the secondary water in a 
UTSG-. The purpose of the introduction of colder feedwater 
directly into this region in the IBtJTSGr design is to increase 
the temperature difference and increase the heat transfer from 
the primary water, 

2.2 Model of the steam generator [6]: 

The steam generator model used in the present study is 
shown in Big. 2.2. The steam generator is represented in two 
sections, preheater (economizer) and an evaporator. In readity 
the steam generators of the have U~tube geometry, the 

present model simplifies it to parallel tube geometry but retains 
all the significant factors such as the total heat transfer 
area, total amount of heat transferred, preheater. Lack of 
detailed data about the internals of the steam generator, such 



Vapour 


Fig. 2 



r.^ Indicate heat transfer 

rT' r,rr,,,t s — Secondary Circuit 

— Primary Circuit „ 

_ Temperature of tube Q - Heat Flux 


Simulation model of Steam Generator 





9 


as the pressure drop due to baffles, quality of steam at the 
outlet of the generator before the steam oseparator etc. is 
circumvented by the pressure balance and heat balance equations 
as described in the later chapters. 

As indicated in the figure the coolant from the reactor 
outlet enters the primary side of the ^nerator as a hot 
stream in the evaporator section and leaves the generator 
from the preheater section as a cold stream to enter the reactor. 
The heat flow from the primary side is conducted throu^ the 
tube walls of the steam generator and enters the secondary 
stream. The figure also indicates that the secondary stream 
leaves the evaporator section of the generator as vapor and 
part of the stream is recirculated. Some make up feedwater 
enters the preheater section and mixes with the secondary stream 
after the preheater. 
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CHiLPTSR 3 

Analytical Ilodel for Simulating the Steam Generator 
3*1 Model for the steady states 


A simplified ma.thematical model and a computer code were 
developed by A.O. Pinto [7] for the steady state analysis of 
a steam generator utilized in a 600 MWe pressurized water 
reactor nuclear power plants. This model takes different flow 
patterns and the corresponding heat transfer regimes into 
account. The details of this model can be found in reference 6. 
However, a brief description of this model is given in the 
foil owing se ctioiis , 

3.1.1 Pundamental equations of conservation; 

In the steady state model [7] it is assumed that the 
flow is vertical. There is no transfer of wcrk across each, 
control volume indicated in Tig. 2.2. The area of flow cross 
section is assumed constant. Thus the three fundamental 
equations of mass, energy and momentum conservation beccme 

Mass = _|| (3.1) 


where W = PAY. 



^ 3 (PAY) 

3 z Tz 


Energy; 


(5.2) 
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Momentum: 


2 

where e = u + ‘V/2 + gZ is the energy per unit 
mass of the fluid. 



A 3iP 

■^ST " Sz 


(3.3) 


3»1j 2 Equations of conservation at stead3^ state: 

At steady state the rate of change of any propertjr with 
respect to time is aero. So the equations '^•1, 3.2 and 3.3 
will become 


¥ = Constant 


3 


.i» 3 


^ dZ 


<L. 

dZ ■" dZ 


(PAY) = 0 


^IZ d( W) 
dZ ^ dZ dZ 


Ap g = 0 


Equation of conservation of energy [7]j 
Using the relations 
Y^ 

e = u + ^ + gZ 

H = u + ~ 

P 

¥ = PAY 


* 

G = ¥/A 


(3.4) 

(3.5) 

(3.6) 


V = 1/p 
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the equatim 3-5 can he transformed to the folloTfing form 
(for details see appendix A): 


dAp 

^ dT" 


- y 


dZ 


¥G^v 


dv 

dZ 


- ¥g = 0 


(3.7) 


Integrating this equation with respect to Z from jxniction n 
to n+1, the limits of control volume n, we get 


Vl 








n+1 


wh €£r e ^ ^ — ^n * 

Control volume n 

heat transfer in control volume n» 


“ ^n^ “ I"* 

0 

is the total amount of 


3.1*4 Equation of conservation of momentum [7]s 

Let us integrate the equation 3*6 with reject to Z 
from the centre of control volume (n-1) to that of control 
volume n. The following two assumptions are made to simplify 
the integration. 

(a) The frictional pressure in one half of any control volume 
is equal to one half of the frictional pressure of the 
whole control volume . 

(b) The velocity of fluid at the outlet of control volume (n~l) 
is equal to the velocity of the fluid at the inlet of 
control volume n. 
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Utilizing these two assumptions and the relationship 
G = U/A, we get (for details see appendix B) the following 
equation. 



1 

2- A 




) 


G 

144- gj 



V. 


n-1 


) 


G 


144 'g. 




(Y^ “ *" l44-g. 


. 4 ^ 


V.. 


(3.9) 


n 


The quantity without har represents the value at the junction 
and the quantity with har represents the average value in a 
control volume • 


An examination of the terms of equation 3.9 will lead 
it to the following familiar form. 


■* 2 frictional^ 


'n 


'n-l 


C.'V.(n-l) 


+ (AP 


fri ctional 


) 


CV(n} 


- toPaooeleratlon>g_^_(^_y ^ “ (''^acoeleratlon>j_^ ^ 

~ ^^^elevation^ (3.10) 

3.1.5 Plow patterns and heat transfer [6]: 

The liquid in the primary heat transport system of the 
steam generator is always in the form of compressed liquid and 
hence the mode of heat transfer is pure convection. On the 
secondary side the liquid enters the generator as saturated 


14 


Id. quid, llie secondary liquid 'Will encounter bubbly flow and 
an-iular flow regimes and the heat transfer regimes will be pure 
convection, local boiling, nucleate boiling, forced convection 
and stable film boiling. The various flow regimes and the heat 
transfer modes and the correlations used in each mode are indi- 
cated in Table 3*1 

3«2 Mathematical Model for Simulating the Transients 

The model developed to simxilafce the transients of the 
steam generator has been described in details in bhe follo?riLng 
sections • 

3 . 2.1 Squations of conservation for analysing the transients! 
The conservation equations are : 


Mass ; 

tl 

- 11 

3Z 


(3.11) 

Energy: 

A 

^ at 

3 Aj 

= q g_ - 

3 (¥e) _ 3 (PAV) 

3 Z 3 Z 

(3.12) 

Momentum: 

A 

3 t 

- - 
“ 3Z 


( 3 . 13 ) 


where e = u + V^/2 + gSis the energy per unit mass of the 
fluid. It can be assumed that the effect of g is negligibly 
small and hence the term gZ in the expression of e and A^g of 
equation 3- 13 can be neglected. 




r~l 

tr\ 

CD 

r-^l 

Bi 
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3.2.2 Equation of conservation of energy: 


as 


A 


Using e = u + V‘^/2 the equation 3.12 can he written 


3(pu+P'V^/2) „ ^ 3(Wu + WV^/2) a (PAV) 

at ~ ^ STg §T” a z 


which on simplification yields 




a V 

+ APV ^ + ^ 


Xi JLI 

2 ’ 3Z 


3Z 


Jel ^ ^ 

(PAV) 


3 (¥u) 
3 Z 


- W 


3Z 


Substituting in ihis equation the value of A ^ we get on 
simpli f i cation 


3(pu} 


Agi 


3t “ ^ az 


ApY 


37 a ( Wu) 3(PAV) 3V 


at 3 z 


3 Z 


•WV (3.12a) 
c-Z 


Erom equation 3*13 we get 


AP + AY H - - A II 


w - v^Ll 
"a-Y " ^3^ -§z^ 


Multiplying both sides by V and substituting the value of 
d o 

A into this equation we get on simplification 

A(.v||= - W||- t’-Se (3.13a) 

Substituting the value at in the equation 3# 12a it can 

be simplified to the following form. 
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A 


a(pu) 

3 t 




3i'- 


E 


3Z 


3 Z 


3Z 


3Z 


(3a2b) 


3.2.? System of equations; 

¥e no-w have the following system of equations to h^ 
used to simulate the transients. 


A _ aSLI. 3^^^) 

3t ~ ”^3 Z " 3Z ~ aZ 


( 3 . 14 ) 


3( P ul 


3t 




3 A^i 

3^“ 


3( Wu) 
3Z 




mgiii -MiMiMxt 

Tz 


( 3 . 15 ) 


( 3 . 16 ) 


Integration of these equations with respect to Z over a control 
volume i hounded by junction i and i+1 will lead to the following 
forms » 


(a) The equation of momentum conservation gives; 


i+1 


i+1 


3 . . , 3(W) 

^ / ApY dz = - / 

i ^ 

^ i+1 

j Ap V dz 


dz 


i+1 

/ 

i 


i+1 

dz ~A / ^ dZ 
i 


Let us define 


dZ 


With this definition we get 
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dV/. 

^ dt 


■^^i+1 - - WiJ - P. 


Ki 


But W = PAY, i.e., W= ¥VpA and 

- AZ 


¥| 


^la ~ % 


2 AP^ P. 
e 1 


So the momentum equation takes the fcrm: 


A = Kt ( i - 


i+1^ A‘ AZ 


1 

-1- ( ^ 


¥? , 

1+1 

Pi+1 


¥7 


% 


2iU), Pi 


(3.17) 


(h) Ihe equation of energy conservation will give: 
i+1 i+1 i+1 

/ Apu dZ = 


X. 

3t 


i+1 

; dZ 

i 


3 Z 


f dZ 

^ 3Z 


i+l 

- / 

i 


1 

i-hl 


9F. 


PA II az + / Y 3^ 


dZ 


I'/e now define 


i+1 

= / A Pu dZ 

i 



i+1 

aA^ 
/ q “^ 2 ^ 


i 
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i+1 

/ 


PA II dZ 


W. ¥. 

11+111 p^ 


and 


i+1 

/ 


d p 

^TT 


^1 ^K,i 


= £ 


8 ¥. ' 

X 


2A^ 

e 1 


AZ 


¥i1;h these simplif 5 ang assumptions we get 


U = 


% + Vi - h+i"i+i 


W-f K- 

. - -i) 

^ Pi+i Pi 


+ f. 




^ 2A^ D P.^ 

e 1 


AZ 


(3.18) 


( c) The equation of mass conservation gives: 


i+1 

3 A ^ 

STt / Ap dz 

i 


i+1 


- / 

i 


3 Z 


dz 


Let us define i+1 

/ Ap dz 




i+1 

A dz 


With this definition we get 


P. ~ A^AZ “ %+l^ 


(3.19) 
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Equations 5»17> 3*18 and 3*19 have to be modified, to incorporate 
the cwo— phase multiplier r and make the units consistent* To 
make the units consistent we have to multiply appropriate term 
by any one of the following three factors. 

2 . 2 

^1 = Sq ^ 3600^ X 144 ^ = 0,600444 x 10^ 

Hr'^ ft*^ 

- 2 

02 = 144 ^ X 0.001285 = 0.18504 

ft 

- 2 

03 = — X ^ ^ X 0,001285 ffe-vv = 0.308172x10“^^ 

^ ^c 3600 ^ Sec^ Ibf-ft 


Eow we have 


in 

P . 

1 


1 

A. A 2 





(3.20) 







+ 


^i ~ \+l ^i+1 


- 1. 





i +1 


+ 


ff AZ 

w JL 


( 3 . 21 ) 





s t;2 


OgA 

AZ 


P ) - 
^ - i+1^ 


f. 

1 


¥4 


2A D, 


— (r.. 

a.az ¥ 







( 3 . 22 ) 
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Now we xjill proceed to express the junction variables in 
terins of control volume variables* 

It is assumed that 


¥. 


(¥. i.e., ¥. 


1 ^ ^“i ^ “i+l^» "i+1 




= i di + i.e., - Uj_ 


^ = i C‘’i + Pi+i). i-a-. Oi+i = 2Pi - l>i 


let us define specific internal energ 7 

i+1 


^i 


/ A P u dz 

i 

i+1 

/ A p dz 
i 



A AZ 


Pressure P. can be expressed as a function of and 


X ♦ 6 ♦ f 


p. = p.C^, P^) 

1 1 i' 1 


These variables will take the following foim at junction 2i 


Wg = 2 Wl - % 

“2 = 2 "l - "l 


P2 = 2Pl - Ol 

Pg = " "i' ^ ^ ~ “b.) 

= P2Cu^, P 2* ^1* ** 1^ ^ 
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At junction 3: 


¥3 = 2 ¥2 - ¥2 


2 ¥2 - (2 ¥^ - ¥^) 


2 ¥2 - 2 ¥^ + ¥3^ 


Similarly, 


= 2 U 2 ~ 2 

p = 2Po — 2P-,+p-) 
^ 2 i J. 


“ P5 (u2» p2» Pi» Pi) 


At junction 4 : 


¥4 = 2 ¥^ - 2 ¥2 + 2 ¥3^ - ¥^ 


= 2 - 2 U 2 + 2 ~ 


^4 - 2 P^ - 2 P2 + 2 p^ - P^ 


^4 “ ^ 1 ' ^ 1 ' ^ 5 ' ^ 2 * ^ 1 ' ^ 1 ^ ^ 1 ^ 


At junction n: 


2 ¥„ , + 2 , - + (-1)“"^ ¥• 


''n = 2 ^ "„_2 ‘ ''n_3 - - v— / »i 


for n = 2, 3» 4, 


, ¥• 


1 .e 


• f 
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\ = (-1)^““^ + 2 I (-1) 


n-1 


n+k-1 




E=1 


K 


for n = 2, 3, 4, 


, F. 


Similarly for n = 2, 3, 4,...., IT, 


rL~l 


= (-1)^"^ ^1 + 2 I (-1) 

E=1 

n-1 

VH -t ^ n— K— 1 

= (-1)^“^ + 2 Z (-1) P 


n+E-1 - 


n 


E 


Ekl 


and 


“ ^n^\L-l» ’^n--2' 


^1» Pn-1»***'» ^1* ^1^ 


Using these expressions the equations 3.20, 3.21 -and 
3.22 can he written as 

n-1 

r{r-i')^’“^ ¥- + 

AA2r 


-■n = h - i, (-I)""""'" 2 \) 


Ei=l 
n 

.n ,, t> , ^ ^n^-E 


-{(-!)"" ¥^+Z (-ir*^^ 2 ly] (3.23) 


k=l 


n-1 

= 4x + [(-1)“"^ '*1 + I (-1)“'"^ 2 ] X [{-1)'^-^ Uj^ 

E=1 


n-1 n 

+ I 2 ;ig.] - [(-1)'^ +J (-1)“+^ 2W^ 

E=1 ‘ E=1 

n 

X [(-1)“ + I (-1)“+^ 2 Ug ] 

E^l 
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n 




E=1 


I (-1) 


n+Z _ B , ^n+E-l 

2 ( i^a-1 „ i 

K (*-1) ¥-l g^-L 


2 ¥n 


E 


n 

I 

E=1 


(-1)^ 1 + I (-1)^'^^ 2 TJg. (-1)^“^ P3_ -2pj^ 

&=1 


,n+E-l, 


“] 


+ 






2A^ D. 
e 



A Z 


(3.24) 


¥ 


n 


On A 

Tz 


^rS ^n~l' %-2» ••*» ^n-l» n-2» V 


AZ 


Wl %~1» \t-2» ••"» Pii» Pn-1»***» ^ 2 .* 

n~l ^ 


+ -r-v^ X 


[(-1)’^“^ + I 2 \] 


E=1 


E-* 


A,AZ ^ n~l 

P, +' I (-1)^+^-^ 2 Pg] 


E^l 


C(-l)^ ¥-1 + E 2 ¥j^ ]^ 


n 


X 


Itrl 


A «/|||^2 ^ -n-i-TT 

[(-!)“ p + E (-1)“+^ 2 Pj. ] 

1 E=1 ^ 


V 


n 


n 


2A D P 
e n 


(3.25) 


5.2*4 Solution of the system of equations: 

Bach control volume is represented hy the equations 3*23 » 
3.24 and 3.25. So for each control volume there will be 3 
equations with 3 unknowns. Porsching [10] had pointed out that 
for any thermal -hydraulic network involving the three basic 
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conservation equations, an explicit numerical scheme can be 
formulated which is consistent and stable for some appropriate 
time-step , 

(a) Derivation of the numerical scheme: 

Equations 4«13» 4«14 and 4*15 can be written as 

* = iP* If 3nd. 3 ~ IjL 

X J 

U = p i = 1, 2^*.., N and 3 = F+1, 21T 

X j 

f . = P- , U = 1, 2,.^,, IT and 3 = 21T4-1, 2F-h2, ♦ • . . ,31T 

X J 

Let us define 

^ = Ve ctor [ ^2 ^ *^2 ^ ^ ^^ 2 * * * * ^ ^IT^ 

and 

p (t^y) = Yector * * * ^^IT^ ^N+l^^IT+2^ * * ^ 

^2IT+l' ^2N+2»*’*»^3 n5 

It can be written using these two definitions 
J = P (t, y) 

So the equation to be solved is 

y = P(t,y) with the initial condition y(o) - yQ^ 



The given differential equation can be witten as 
dy 

df = I y) 

i.e., At 

= ■f ^t 

= { + A } A t 


If we approximate 



dP 

«ln " 

Ay^^"- 


then 


dt ^ 


So we 

now have 




,,h+l 

A I 

[ p“ + If 1 Ay^"^^ 

] A t 

X « 0 • jp 

[I - ^1 
^ ay 1 n 

At ] Ay^-"^ = f 

At 


The matrix representation of this equation will be 
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~A t ^ 

3tj. 

- At ^ 

a 


ail 

l-At — ^ 
aUi 

aC 

- At *—3^ 

at. 

— At 

3Pi 

a t. 

— ^ t 

3 0, 

a ¥. 

i-At — ^ 

a 


t 

1 


! ^ . 

1 i 

! 


1 

# 

I - 

P- 

1 



! 

f 



: 

A\ 

f 

i 

] 

♦ 

h 

( 



J 


A t 


Before evaluating these derivatives let us first estab- 
lish a few identities which will be used for the evaluation 
of the derivatives. It has been shown that 


^n - 


n-l 


n-1 

+ 2 I (- 1 ) 

Erl 


n+K-1 


K 


ac 


and 


3Pn 
9 Pj 


aimilarly, 




(^1)^+3-’! .2 


(-1)^+3 , 2 


and 


3U. 


8 u 


3U;. 




i-l 


^ = (-l)^'*'^' • 2 


/ 
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ilso 


U. 


X .e. , 


u . 


3 u. 




A 


3 U;j A p|j A z 


Using tliese identities and talking the derivatives of equations 
3 . 23 , 3*24 and 3»25 with respect to P . , U. and ¥. for n = i 

ill 3. 3* 

and utilizing the following three notations. 


m-1 


m-l 


= 

m 


(-1) ¥-, + 2 I (-1) 

E=1 


m+K-l 


¥ 


K 


m-1 


m-1 


u 


m 


(— 1) Un +2 I (“*1) 

I &=1 


m+E-1 




= 

m 


m-1 

(- 1 ) *^1 + 2 


m-1 

I 

£=1 


(- 1 ) 


m+E-1 


Re 


¥e get 


3p . 


= 0 


3pj 

3U. 


= 0 


3 ¥ 


A-AZ 
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9 









3?- , W. .. ¥, 

~ Op ( ~ ) 


Pi 


Pi+1 


+ 02* P * 2 ' 


\+l 

"i+1 


- 


G 


-3 

¥. 

1 

X 


az 





A-Az* 


p 

1 


c. 


(^>1 ( 

3u ^i+l 


)- 



) 


2u 


i+1 


2 P. 



1 

2 


h t az 

A^*D ‘P.^ 

e 1 



J- 

AZ 




-'i+l 

A*A2 


¥^ 

^i+1 


+ 




2 ADe 




_2 ( ^ 

AZ^ Pi 


f . ¥i 


^i+1 ^i+1 


A'D 


e 


Pi 


A» A Z* 


Pi+1 


¥sing these values in the above 3x5 matrix the values of 

^ ATJ> and a¥- can be calculated for each control volume 
1 i ^ n-l -n+l -n+l 

and knowing these values f 0"^ and ¥i can be 


determined by using 
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_ n+1 

P = P + 

i i i 


n+l 



n 

U. + A U. 

-L. JL 


_n+l _n 
% = \ + A 


-n+l _n+l _n.+l 

Once P _ » U and W are known, oth-cr paxameters 

i i i 

can be easily determined* 
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CHAPISR 4 

DBS GRIP TIOrr OP THE OOMPUIM GODS 

0 

4.1 Input Variables: 

The input to the computer code is fairly simple. Some 
of the data can be directly taken from the steam table. But 
other data have to be generated by polynomial regression using 
the data of compressed liquid (water) . 

The required input data cards are listed below. The units 
of each parameter and their input format are also included. 


Card Fo. Variable 

Value 

Format 

1-118 AGUiS (1,1) 

Saturation pressure 
in psi 

P 10.4 

AGUAS (1,2) 

Saturation tonperature 
in op 

P 10,4 

AGUAS (1,3) 

Sp. Vol. of liquid in 
ft ^/ibm ♦ 

P 10.4 

AGUAS (1,4) 

Sp. Vol. of vapour in 
ft ^/Ibm . 

P 10,4 

AGUAS (1,5) 

Sp. enthalpy of liquid 
in Btu/lbm 

P 10.4 

AGUAS (I»6) 

Sp. enthalpy of vaipour 
in Btu/lbm 

P 10.4 


The cards are in decreasing order of saturation pressure. 
119-126 ((AGUAG (I,J), J = 1, 12), I = 1, 4) 
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Pormat is 6S12,6 for each card and there 
■will he 6 values on each card. 

J = 12 different values of temperature 
1=1 Temperature 
1=2 Constant part 
1=5 Coefficient of P 

1=4 Coefficient of 

127 PRBGrB Limits of quality regions 6P 10.4 

6 values 

128-166 C0BPB(I,J) Interpolated values 5E 15 *6 

((C0BPB(I,J), J = 1, 15), I = 1, 13) 

5 values on each card. 

1=15 Quality for the curve 

J=l,15 Indicates the pressure 

set to be used. 

1=1-6 A for each set 

1=7-12 B for each set 

(G0EPB(15,E) , &=1,15) Quality of each 

of iiie 15 curves 
of Baroczy plot. 

167-168 G-AMAT Gamma factor 7P 10.4 

((GAMAT(I,J) , J=l,7), 1=1,2) 

J=1 


J=2 


Indicates pressure 
Indicates gama factor 
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Card ITo . 

Variable 

Value 

Format 

169 

GNOME 

Name of th.e programme 

20A 

170 

NVG 

No* of control volume 

13 


NiyiGL 

Maximum No. of 
iterations 

13 


lESGR 

Output option 

= 0, Print final value 

=1, Print all iterations 

=2, Print all details 

11 


Pi®CQ 

Precision required 

PIO .5 


GRAV 

Acceleration due to 
gravity 

PIO .5 

171 

NTP 

No. of TJ-tubes in primary 

17 


MALHA 

Type of pitch 

= 1, Square pitch 

=2, Triangular pitch 

11 


ATTO 

Total heat transfer area 
(ft^) 

P10.2 


DBZ 

External diameter (inches)P10.2 


DIN 

Internal diameter (inches) 

P10.2 


PASSO 

Spacing between tubes 
(inches) 

P10,2 


COBPP 

Perm factor 

P10.2 

172 

PPJ(l) 

Inlet steam pressure 
(Psi) 

P10.2 


TPJ(l) 

Inlet steam temperature 
(°F) 

P10*2 
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Card Ho. 

Variable 

Value 

Format 


HPJ (1) 

Inlet steam enthalpy 
(Btu/lhm) 

F10.2 


PSP 

Estimated primary 
pressure (Psi) 

F10.2 


w 

Inlet flow rate in 
primary (Ihm/hr) 

B12.6 

173 

PEM 

Pressure of feed water 
(Psi) 

F10.2 


IBM 

Temperature of feed 
water (°E) 

F10.2 


HEM 

Enthalpy of feed water 
(Btu/lbm) 

F10.2 


PSS 

Exit pressure (Psi) 

F10.2 


XExit 

Quality at exit 

F10.2 


HS 

Plow rate of feed water 
(Ihm/hr) 

B12.6 

174 

QTOT 

Total heat transfer rate 
(Btu/hr) 

B12.6 

175 

OP 

Calculation option 

= 1, calculate only the 
steady state values 

/ 1, calculate both 

steady state values 
and transients. 

11 

176-183 

OP7AGCl,J) 

Pol 3 aiomials expressing 
internal energy (u) as a 
function of sp.vol.(v) 

5X, F7.1, 
5Z, 3B16*5 


((UJVAG(I,J) , cr=1^4), 1=1,8) 
4 values on eadb. card. 
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Card No. Variable Value 


Format 


184-191 


192 


1=8 different values at 

pressure 


H 

II 

Pressure 


J = 2 

Constant pai-t 


it 

Coefficient of v 


II 

Coefficient of v^ 


The cards are 

in increasing order of 


the pressure. 

These are for compressed water. 

VHJAG(I,J) 

Polynomials e:^ressing 
sp.vol,(v) as a function 
of internal energy (u) 

5X, 17*1, 
5Z, 3E16.: 

((VJUAG(I,J), 

J=1,4),I=1,8) 


4 values on each card. 


1 = 8 

different pressure 
values 


J = 1 

Pressure 


J = 2 

Constant part 


11 

Coefficient of u 


J = 4 

2 

Coefficient of u 


The cards are 

in increasing order of 


the pressure. 

These are for compressed water. 

Time (Sec) 

]}uration over which the 
transient has to be 
cailculated. 

No format 
is needed 

NIT 

No. of intervals 

No format 
is needed 
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Card No. Variable Value Pormat 


193-207 VPPH (l,J) Polynomials expressing 3X, P6.1, 

sp.vol. (v) as a 3S19.7 

function of pressure (P) 
at constant enthalpy (h) 

((VPPH J=l,4), 1=1,15) 

4 values on each card. 



I = 

15 

different enthalpy 
values • 




J = 

1 

Enthalpy 




J = 

2 

Constant part 




J = 

3 

Coefficient of P 




J = 

4 

Coefficient of P^ 




The 

cards are 

in increasing order of 




the 

enthalpy. 

These are for compressed 

water . 

208 

VBPPJ (2) 

S.p.Vol, of compressed 

No 

format 




water at the inlet 

is 

needed 




condition of the pre- 
heater of primary 
circuit • 



209 

VESPJ (1) 

Sp. volume of water at 

No 

format 




the inlet condition of 
the preheater of second- 
ary circuit. 

is 

needed 


4.2 Steady state calculation: 

Por calculating the steady state values, i.e., the starting 
■values for transient calculation the conputer code and the 
model developed by A.O. Pinto [7] are used. 
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4.2.1 Plow- chart i 

The procedure involved in the above computer code is 
given in the form of a flow-chart in figure 4.1. The function 
of various subroutines will be discussed in the next section. 

4.2.2 Subroutines [7]: 

(a) QPAG-C - It calculates specific heat of compressed 

water at constant pressure as a fmction of temperature (°P) 
at various pressure range. The unit of specific heat is 
Btu/lb-^P. 

(i) Op = 0.14931.10“^ T^ - 0.45187.10“^ T + 1.0241 for 

P < 1500 Psi. 

(ii) 1500 Psi : P C 2500 Psi 

= 0.90724.10“^ T^ - 0.23174.10“^ T + 1.0049, T^ 300°P 

Gp = 0.38132.10“'^ T^ - 0.453.10“^ T^ - 0.18281.10“^ T - 1.4375 
for T > 300°P. 

(iii) 2500 Psi < P ^ 3500 Psi 

Cp = 0.90941.10“^ T^ - 0.24668x10“^ T -h 1,0028 for T 300°P 
Gp = 0.47072,10“^ T^ - 0,31944.10“^ T + 1.5651 
for 300 < T ^ 620 °F 

Op * 0.22404.10“^ T^ - 0.10377.10“^ T^ - 0.13194 T + 69.772 
for T > 620 °F. 
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(l)) GPAGS — It calculates specific heat Cp of saturated 
water as a function of temperature T* 

Op = 0.10313.10"^ - 0 , 25129 . 10 "^ T + 1.0125 for T 410 °P 

Op = 0.13446.10"^ 0.11339.10"^ T + 3.4878 

for 410 ^ r <■ 630 °P 

Cp = 0.21005 !D - 11.543, 630 ^ T < 660 °P 

Cp = 0,58005 T - 35.97, 660 ^ T < 680 °P 

Op 0.20068 T - 132.97, T '> 680 °P 

( c) CTAGrC - It calculates the conductivity of compressed 
water in Btu/hr-ft-°p as a function of press \ire (psi) 
and temperature (*^P) . 

(i) P£. 1500 Psi 

k = 0,1131.10"^ + 0 . 64326 , 10 "^ I + 0.30963 

(ii) 1500 < P :< 2500 Psi 

k = -0.11059.10“^ + 0 . 63854 . 10 "^ T + 0.31182 

(iii) 2500 C P £ 3500 Psi 

k = -0.114 95 . 10"5 + 0 . 6 7117. 10 "^ T+ 0 . 3105 . 

(d) CTAGrC -• It calculates the conductivity (k) of saturated 

water as a function of temperature (T). 



k = -0.24 17 3. 10“^ + 0.99835.10"’^ T + 0.28862 


for T £ 100 

k = -0.10408.10*"^ + 0.5628.10"*^ + 0.32132 

for 100 < T ^ 400 °E, 

k = -0. 20081. 10""^ + 0 . 15749 . 10 "*^ T + 0.066 24 

for T > 400 °E 

OTMTP - It calculates the conductivity (k) of the material 
of the tube, Ihe material chosen, in this problem 
is Inconel-718. 

k = 57.78 (-1.98x10"*^ + 8.91x10“^ T + 0.151) 

where k is in Btu/hr-ft-°E and I in °G. 

PTMDP - It calculates the multiplication factor needed 
to evaluate two-phase frictional pressure drop. 

Eor this Baroczy's plot [9] is used. The Z-axis 
of this log-log plot is divided into 6 regions 
and the segments in each region is assumed to be 
linear such that 

0^ = Bxp (A + B In I). 

Lo 

A and B are determined using the plot. The 
SIX regions arei 
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legion 

I ; 

0.00001 <- 

I 


0.001 

Region 

II : 

0.001 < 

I 

'C 

0.003 

Region 

III : 

0.003 < 

I 

c 

0.007 

Region 

IV ; 

0.007 < 

I 


0.03 

Region 

Y ; 

0.03 -- 

I 


0.1 

Region 

VI : 

0.1 .s-l 

I 


1.0 

where 

I = ( 

(- 

-L V" 



( QPOCV — It calculates ]?, S and o to be used m Chen’s 
correlation. 


F = B^qp 

[0.60358 + 0.21681 ln(X~J)], 

0.1 


0.2 

F = Bxp 

[0.94034 + 0.42604 ln(Z“^)], 

0.2 

^t 

0.6 

F = Bxp 

[1.0359 + 0.61316 In (X~^)], 

0.6 

^-1 

^t 

2.0 

F = B2cp 

[0.96332 + 0.71791 ln(X~^)], 

2.0 

..-1 

^“•tt 

100,0 


S = 0.433 21. 10~^° - 0.10254.10"'^ B + 0.97877 

For 14990 

S = 0. 54595. 10~^^ B^ - 0,33653.10**^ B + 0.65916 
For 10^ B f 3.0x10^ 

S = 0.10265.10"*^^ B^ - 0.99733.10*“^ B + 0.34334 
For 3x10^ c B ^ 5x10^ 
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a = 0.11075.10”^^ - 0.18358.10”^ B + O.I64O6 

For 5x10^ .i: B < 7x10^ 

S = 0.9, For B j> 7x10^ 

where B = Re^. 

0 = 0.157.10”^ - 0.36872.10”^ P + 0.0033934 

Per P < 1000 Psi 

0 = 0.2023.10“^ P^ - 0,1402.10”^ P + 0.0024308 
Por 1000 -< P £. 3000 Psi 

where 0 is in Ibf/ft. 

(h) TLSUB - It determines the temper atnre of compressed 

water for a given pressure ( P, psi) and specific 
enthalpy (H, Btu/lbm) . 

H = 0.26127.10”^ P + 68.0349, T = 100 °P 

H = 0.22714.10”^ P + 167.946, T = 200 °F 

H = 0.18834.10“^ P + 269.445, T = 300 °P 

H = 0.49075.10”'^ P^ + 0.10506.10”^ P + 374.7088, T = 4OO °F 

H = 0.10528,10“^ P^ -- 0.56072.10"^ P + 488.1221, I = 500 °P 
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H = 0.81068,10 ^ - 0 . 84505 . 10 "“^ P + 628,158, 1 = 600 

H = 0.13558.10'*^ P^ - 0,13685.10”^ p + 666.845, T = 620 

H = 0.18628.10”^ P^ - 0.20086.10”^ P + 711.776, T = 640 

H = 0.53821.10”^ P^ - 0.46992.10“^ P + 794.942, T = 660 

H = 0 . 75397 . 10 ”^ P^ - 0.7298,10“^ P + 896.170, T = 680 

H = 0.19416.10”'^ p2 ~ 0.16793.1* + 1104.3, T = 690 

H = 0.24845.10”'^ P^ - 0.22748.P + 1274.58, T = 700 °P. 

) VIAGO ~ It determines tlie Yiscosity (ji) of compressed 

water in Ibm/Hr-ft as a function of pressure (p, 
psi) and temperature (T, °p) . 

(i) Prf 1500 Psi 

p = 0.3705.10”^ - 0.87236.10”^ T + 6.6531 

for T < 110 °P 

\x = - 0 . 66532 . 10 “'^ lP + 0,66982.10*"^ - 0.23742.10”^ T 

+ 3.3519 for 110 ^ 400 

p = 0.14575.10”^ - 0,20341.10“^ I + 0.90965 

for T >■ 400 
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(li) 1500 <CP < 2500 Psi 

[i = 0,36^96.10 ^ - 0,86013*10“^ 1 + 6.5907, T ^110 °P 

= -0,65878,10 T^h- 0,66502,10“'^ 0,23542.10“^ T 

+ 3*3435, 110 °P < T < 400 °P 

[I = 0,11304*10“^ - 0.17363.10“^ 1 + 0.84679, T>400 °P 

(lii) 2500 Psi < P <: 3500 Psi 

[X = 0.35894.10“^ - 0*84719.10“^ T + 6.5258, T < 110°P 

tx = -0.65076.10”'^ lP - 0,65502,10“^ T ^- 0.23311*10“^T 

+ 3*3325, 110 < T < 400 °P 

= 0.52496,10”^ - 0,11349*10”^ T + 0,70403, T >400 °P 

(0) VIAGS - It determines the viscosity of saturated water 

([j.) in Ib/lir-ft as a function of temperabure ( 1) * 

|i = 0.37664*10”^ - 0.88541.10“^ T + 6.7181, I ^110 °P 

|x = -0*67639*10”'^ + 0*68025*10”'^ - 0.24029.10”^ T 

+ 3.3665, 110°P T < 400°P 

\i = 0.29905*10”^ ip- - 0.87608.10”^ T + 0.62286, I > 400°P 

(k) VIViS - It calculates the viscosity of satxxrated vapour 
in Ibm/Hf-ft as a function of pressure (psi) . 

p = 1. 15826 (-0.243 53 aO”® + 0*13778.10”'^ P + 0.029109) 
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4.3 Transient calculation: 

The details of the oomputei* programme deyeloped for 
simulatjjig the transients of the steam generator has been 
discussed in the next section. 

4 . 3.1 Plow chart ; 

The logic of the progi^am is presented in the form of 
a flow diagram in figure A, 2, 


4 , 3.2 Subroutines: 

(a) SOLUC - This subroutine calls a few other subroutines 

and simulates "Ihe transients . The function of 
this subroutine is given in the form of a flow 
diagram in figure 4*3 • 

(b) ICRIlIl - It determines for both compressed and 

saturated water. To evaluate this the following 
relationship is made use of. 


{ ^ ) 


u 


1 


>v 'u 


(i) Compressed water: 

For compressed water the specific volume (v) was expressed 
as a function of internal energy (u) at different pressures 
using the data from the compressed water table of fief, 11. The 
polynomials are of the form 
V = a + bu + ou 


where a, b and c are: 



45 


Pre ss ure 
(Psi) 

a 

b 

c 

400.0 

0.15991E-01 

O. 902 IIE -06 

O.I 64 O 8 E -07 

600.0 

0.15979B-01 

0.92847B-06 

0.16303B-07 

800.0 

0.15987E-01 

0.40455B-06 

0.17907B-07 

1000.0 

0.15974E-01 

0.46761B-06 

0.1769aB-07 

1500.0 

0,15945E-01 

0.58317E-06 

O.I 7275 E-O 7 

2000.0 

0.1604^B~01 

-0.20277E-05 

0,23350B-07 

2500.0 

0.16176B-O1 

-0.44908E~05 

0.27727B-07 

5000.0 

0.165O3E-01 

-O. 953 O 8 B-O 5 

0.35265B-07 


To evaluate tlie required derivative the range of pressure 
between which the known pressure exists is first determined. 
Lot us assume that the lower limits is and upper limit is 
1“^* Gone* e spending to these two pressures the specific volumes 
''^lo "“"u determined using the appropriate polynomial* 

The required derivative will be 


( ) 


^10 “ \ 


(ii) Saturated water: 


Por saturated water two other pressures are selected 
such that they are on either side of the known pressure, 
equidistant from the known pressure and do not differ very 
much from the known pressure. At these two pressures the 
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specific volumes and. lienee flie densities ane determined and 
blien the required derivative is calculated by eicpressing it 
in finite difference fona. 

The error introduced by this technique is partly due to 
the fact that the slope of the chord differs sli^tly from 
the slope of the tangent at the mid-point of the arc, and 
partly due to the small difference between adjacent variables. 


( c) DPRjaU - 


This subroutine determines ('^)^ • The tedanique 
followed in this subroutine is same as that des- 
cribed for the subroutine DPIiBR. But in this 
case the internal energy (u) is expressed as a 
function of specific volume (v) for the compressed 
water. The polynomials are of the form 


u = a + bv + c V 


where a, b and c are : 


Pressure 

(Psi) 

a 

b 

c 

400.0 

- 0.14297B+05 

0.154993+07 

- 0.409513+08 

600.0 

- 0.14553B+05 

0.155753+07 

- 0.411793+08 

800.0 

- 0.82765B+04 

0.848103+06 

- 0.205543+08 

1000.0 

- 0.85259E+04 

0.854213+06 

- 0.207163+08 

1500.0 

- 0.844653+04 

0.869153+06 

- 0,211613+08 

2000.0 

- 0.471793+04 

0.450913+06 

- 0.951553+07 


2500.0 

3000.0 

(d) 


“ 0.34271E+04 
“• 0.25480i;+04 


O.3170OB+O6 

0.22901il+06 


“ O.6II8OB+O7 

- 0.3981213+07 


This subroutine determines the quality of steam 
using ITewton-Rapshon’s method. It utilizes the 
following four polynomials which are obtained 
hy fitting curves for the saturation pressure 
of water. The following polynomials are valid 
in the range 401. 7 < p 2532.4 psi. 


v^ _ 0.1745921x10 + 0,5229798x10’"^ P - 0,1625896x10"'^ 

+ 0.5351606x 10"’^2 p3^ 


Vg = 0.1S67797 x 10-"1 - 0.?7 10660x10-2 I + 0.1436388x10-5 p2 
~ 0.2599688x10“^ P^. 

I'l 

Up = 0,3178337x10^ + 0.3054259 P - 0.1020149x10“^ p2 
+ 0.17564660x10“'^ P^, 


u^ = 0.1120552x10'^ + 0.4207413x10“^ P ~ 0.1494585x10“^ P^ 
- 0.4154948x10“^ P^ 

The quality Z = -^—1 = 

\ ~ ^f ^g “ '^f 
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In the Ii.H.S, of th-q ^ 

1-n s e..pre.Sion u and t ai-e too™ and hanoe 

GOn„tant. All variablas other than these two are function of 
presonre. So pressure can he calculated fro. this equation. 
Por this purpose JTewton-hapshon. s method xs ueed. Once 

P P© S S 1X2?© 1 3 IQI O TlaTH D n T 

°an he calculated from the above 

polynomials and hotice thp nnoTu -(-tt- ~r n 

u ucnce me quality L can be determined. 




ARRANGE VALUES IN TAB 


ULAR FORM 


PRINT GENERAL DATA , CA L CULATi 
GEOMETRICAL PARAMETERS 



DETERMINE THE AXIAL D15TPIB 101 
OF TEMPERATURE AT-jD QUALIT-' OF 
VAPOUR IN THE EVAPORATOR 



COE^Rf t 

T DF 

(vEOMET 

IHCAI 

AND 


OE't RAI I 

T-G 

PARA ME 

TE R'; 


CALCULATE THE POWER GENERATED 
IN THE STEAM GENERATOR AND 
COMPARE WITH INPUT POV/EF^ 
RECEIVED FROM REACTOR 


^ IS THE ^ 

REQUIRED PRECISION 
OBTAINED^i.--^ 


DF Tt RMINF THE 
AXIAL D,S‘'[UR JTION 
or PRE S5UPE 8, 

Tf MK IN fH<lMARY 

secondary 


HAS IT 

EXCEEDED MAX.NUMBER OF 
ITER AT ION 2..—^ 


YES 


CALCUIATE THE PARAMETERS OF THE 
PREHEATER. HEIGHT, HEAT FLUX, 
PRESSURE, TEMPERATURE ETC- 


Fig*4'1 Flow chart of the computer programme 
calculating the steady state values; 
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CUi.Tr^OL VO^J.KfS AfiL A-^ 'aLI 

J!'iAC;!orr> a"" i o 


iNiriALlZh THE COf.TRCL VOLUME 
ONF-OR TWO-PnASE 
UTARTh Ar.D ENDS 


INiClALiZE THE T,me AhD CAuCuLA 
THE TIME STEP 


DO 10 u : , r;,T 
NiT= TOTAL NO OF STEPS 


01 T E R M ; N E W h !C h CON T RO ^ VO LU Ml h A S 
ONE -OR TWO- PHASE FLOW AT TIME 


(ALL QMODC ANDCALLULATE "'Ht 
hFAT FloX 


DE^eRM LETHE AMOjNT OF hEA"" i 

^RAfiSFr RRif.G IN EACr CONTPOu VO I JME I 


DLTFRMif.F THE AMOUNT OF HEAT 
TKANbl E RR NO IN THE PREhEATEP 


DETERMINE THE MODE" OF F LOW ir. '^HE SECONDARY CiRCUlT 


EVALUATE w, /5, u , H 0,P, / AT THE INLET OF P^^IMARY LOOP 


CALL SOLUC AND CALCULATE THE VARIABLES LN ALL CONTROL 
VOLJMFS AND AT ALL JUNCTIONS AT TIME FOR PRIMARY 
OF • EVAF-'OPATOR 












SCLUC 

CALCULATION FOR PRIM>^RY OF 
PRt heater 


JSF COf.SERVATIOr; OF MASE & 
ENERGY TO CALCULATE ^FE 
VARIABLES AT THE l^4LET AT 
SECONDARY CiPCUtTSOF PRE- 
HEATER AND EVAPORATER 


SOLUC 

CALCULATION FOR THE 
SECONDARY OF PREHEATER 


SOLUC 

CALCULATION FOR THE 
SECONDARY OF EVAPORATOR 


STORE THESE VALUES 


0 CONTINUE 



PRINT THE RESULT 


'f * ' 















DP REU 

CALCULATE (lP)o IN 
9 u I 

ALl. CONTROL VOLUMES 


DETERMINE THE JUNCTION 
AT WHICH ONE-OR TWO PHASE 
FLOW STARTS & TERMH'.ATES 


D P REU 

CALCULATE ^)pAT ALL JUNCTIONS 


D P R E R 

CALCULATE (|5)n IN ALL CONTROL 
VOLUMES 


D P R E R 


CALCULATE 


)^J ALL JUNCTIONS 


THERE SINGLE PHASE 
FLOW? 


m 


F F R U F 

CALCULATE THE FRICTION FACTOR 
FOR TWO PHASE FLOW 


THERE TWO -PHASE 
FLOW? 


F F R D F 

CALCULATE THE PRESSURE DROP 
DUE TO FRICTION 













THERE ANNULAR 


SET THE 
MULTIPLi GATOR 
R = 1 AT ALL 
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R T R A N 

CALCULATE THE MULTIPLI CATOR 
R AT ALL JUNCTIONS 


R E SOL 
SOLVE FOR 
Ap,AU AND 
Aw S< CA LC- 
ULATE PU ^ 
w AT 


MOr.T G 

CALCULATE THE COEFFICIENTS OF 
MATRIX 


CALCULATE u AND v AND 
DETERMINE /^w,u AT 
ALL JUNCTIONS AT 


P E R P T 

DETERMINE PJ,H AND X IN 
ALL CONTROL VOLUMES 


NEWTON 
CALCULATE P 
AND X 










■w 

M 

■ . -fi -v/ 


DETEMINE P,T,H AND X 
AT ALL JUNCTIONS 


RETURN 


END 


Flow chart of the subroutine 50LUC 


If '.''i ^ ' 

>?■«<{< J ^ r I T- ^ H. 
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GHiPTSH 5 
RESULTS 

5*1 Results and Discussions: 

This chapter deals with the results obtained using the 
computer code developed in this present work. The results will 
be discussed in details in the following two sections. 

5*1.1 Steady state values: 

The steady state values whj-ch are required for transient 
analysis were obtained using the computer code developed by 
A.C. Pinto [7]. These steady state operating values which 
ibolize the steam generator chosen for transient analysis in 
the present work are given in appendix C. 

5*1.2 Transient values: 

The results of the transient analysis are represented in 
graphical forms. Pigures 5*la throu^ 5»6d indicate the tran- 
sient behaviour of the steam generator under perturbed conditions. 

Figures 5.1a through 5.1d represent the change of varia- 
bles at the exist of the evaporator of iiie secondary circuit due 
to a 5,0 X increment in the mass flow rate at the inlet of the 
primary circuit. This increase in mass flow rate in the primary 
circxiit will enhance the heat transfer from the primary to the 
secondary circuit. So the temperature in the secondary circuit 
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will increase and this will raise the q[U.ality. Increase in 
quality means that there is more vapour and 1iie mass flow rate 
of vapour is higher. Thus the pressure will rise because there 
is more vapor and the total mass (liquid and vapor) flow rate 
tCLll in ere ase . 

Tigures 5.2a through 5. 2d indicate the effect of reducing 
pressure hy 30.0 at the inlet of the primary circuit. When 

pressure is reduced hy such a large value there will he two- 
phase flow in the primary circuit and this will increase the 
heat transfer coefficient from the hulk of the fluid to the 
tube wall. This essentially will lead to an increase in heat 
transfer to the secondary circuit and thus the same phenomena 
as in figures 5.1a through 5. Id will occur. 

In figures 5.3a through 5.3d the effect of increasing the 
internal energy of the primary liquid hy 5.0 % has been shoxm. 
Since the internal energy of the primary liquid is high, its 
heat content will also he high and thus the heat transfer from 
the primary to the secondary circuit will increase. So a 
similar effect like figures 5.1a through 5.2d can he observed. 

Figures 5.4a through 5.4d represent the effect of reducing 
feed water flow rate by 5 at the inlet of the secondary 
circuit. Since the flow rate is reduced in the secondary loop 
the amount of water receiving heat from the primary circuit is 
less. So the water in the secondary circuit will have higher 
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temperature and the quality will increase. Thus there will he an 
increase in pressure and flow rate. But the disturbance is very 
small » only 5*0 / » So these changes at the variables at the 
exit of the evaporator of the secondary circuit will also be very 
small. 

ligures 5»5a through 5.5d indicate the effect of reducing 
pressure by 30.0 /C at the inlet of the secondary circuit. The 
reduction in pressure will increase the quality and hence the 
mode of flow of the fluid in the secondary loop will be such 
that the heat transfer coefficient from the tube wall to the 
bulk of the fluid will increase which means an increased heat 
transfer rate. So the temperature will rise which in turn will 
affect the quality and this inciease in quality will increase 
the pressure and mass flow rate. 

ilgures 5*6a throu^ 5»6d indicate the effect of increasing 
the internal energy of the secondary liquid by 5«0 JC. This 
increases the total heat content of the secondary fluid and 
hence the same phenomena as in figures 5 •5a through 5«5d can be 
observed. 

5.2 Conclusi ons i 

In the present study it has been found that a 5«0 % change 
in the feed water flow rate ab the entrance of the secondary 
circuit of the evaporator does not affect the ftanction of the 
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steam generator very much. This trend wa^ established by 
A.C. Pinto [7] in his attempt to find out how the steady state 
values at the exit of evaporator of the secondary circuit vary 
due to different feed water flow rates. 

In all other cases the change in the variables at the out- 
let of the evaporator of thesse condary circuit for some distur- 
bances at the inlet conditions la large. This signifies that 
the mass flow rate, pressure and internal energy at the inlet 
of the primary circuit and the pressure and internal energy 
at the entrance of the secondary circuit play an important role 
in the operation of the steam ^nerator. 

If on© ©xarain®® the graph® csrefnUy one win fin* tliaft 
there are fiiaoontlnnities in ^vost all the eraffee at 0*4 aeo* 
fhl® he due to Bcraerlcal inetabillty* f&# ti*# atep 
th:r the pres®it worls i® 0*001 see. 9© farther ^©011011:10^ 
he SBt&B ifith tiis® step th^iih thla win take imt# 

of GOTmiitmr tine* 



Pressure(psia) 


Pressure at the exit of evaporator 
of secondary circuit 

2800 r 



Time(sec) 


Fig. 5-1a Mass flow rate at the inlet of primary 
circuit was increased by 5-0“/o 




Temperature ( F,' 


Temperature at the exit ot evaporator 
of the secondary circuit 

'"(TV 


I 

SriOl- 



Fig. 5 1b Ma; ; '.;)'V re e at the inlet of primary 
rtrr ins creased bv 5 0% 



lube (Ibm/hr) 


MOSS flow rote Qf the exit of evaporator 

of the secondary circuit 
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•Id Moss flow rots of th© inlet of the primary circuit was 
increased by 5-0°/» 





Quality at the exit of evaporator 
of the secondary circuit 
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5 *2(1 Pressure was reduced by 30-0% at the inlet of 
evaporator of the primary circuit 
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Mass flow rate at the exit of 
evaporator of the secondary circuit 
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il'.Rg. 5'3c The internal energy of primary liquid 
stf'-' was increased by 5-0°/o at the entrance 


of thp evaporator 
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the secondary circuit 

72l ; 

j 

1 

I 

! 

680r 



Time (sec ) ' 


Fig.5-4b Mass flow rate at the inlet of the 

ndary circuit wc- -'xiuceci by 5 07 
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5'5a Ti e pressure at the inlet of the 

S '::Dndary circuit was reduced by 30-0°/< 
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iPPjilTDIZ A 

¥e have 

_ 

e - u + ^ + g2 

2 = u + p/p 

¥ = PaV 

a = ¥/a 

V = 1/p 

The last two terms of eaDa+inn x c: i 

ux equation :5,5 can be written as 

<iC¥e) , d / . 

IZ dZ 


= (ffi) 

= + (?)} 


ITow 


e = u + ^ + g2 = H 


P 

P 


7 

+ ^ + gZ 


(a) 


* ^ 

* . *dZ 


dZ 


dZ ^0^ + ^ 


dZ + s 


de 

dZ 


^ (S 


dZ ^p' 


^ + g + V 32 
dZ ^ ^ dZ 


(b) 


But 7 U = JL ^ fiL ) _ 
dZ ^ dZ '^PA' ~ 


= V 


2 ^ 
dZ 


I (If ^ fl) 

P dZ 


(c) 
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Prom equation (a), (t) and (o) 


we get, 


d(¥e) 

^ (PA7) = V (41 + ^ ^ n2 dv 


^ + g + V e- ) 


So equation 3.5 will become 


.. 


dH 


dZ ^ (12 “* ¥ G-' 


2 ^ 
dZ 


V rr- - -ir/g _ Q 
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^PSITDIZ B 


Equation 3.6 is 


- A ^ - 
dZ 


dp. 


K 


dZ 


d(w) 
dZ 


- Apg = 0 


Eow let us intepnate t>n <5 4 . 

Ills equation with respect to Z from 

the centre of control volume (r\ i ^ +« -f-i, 

(n 1 ) to the centre of control 

volume n. 

The first term will give 


centre at 
C.V.(n) 


/ 


dZ 


dZ 


centre of 
C.V.(n-l) 


(P) - (P) 

centre of centre of 

C*V.(n) C«V«(n—l) 



where p - (p) 

^ at the centre of 

C.Y.(n) 

The second term will become 


J.n. 


djg + 


centre of 
C.7.(n-1) 


eenure of 


C.V.(n) 


dP. 


dZ 


K 


dZ 


centre of 
C.T.(n-l) 


centre of 
G.Y.(n) 

/ 


J .n 
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1 

2 


[ 1 


E,n-1 


+ 



] 


'J.n ’ repx-GSenfcs a unction n. 


Prom the third term we get 


centre of 
C.V.(n) 


/ iLM 

dZ 

centre of 
C.V,(n-l) 


dZ 


J .n 


^ f dV 


+ W 


centre of 
C.7.(n~l) 


centre of 
0.7,(n) 


/ dV 


J .n 


¥,(Y^ 


3’'t J •!! at the centre 
of C.Y.(n-l) 


The last term will be 


+ ¥.(? 


at the centre 
of C.T.(n) 


- Y ^ 

at J.n 


come 


centre of 
O.Y.(n) 


Apg dZ = lpgAZ= 

V 

centre of 
C.Y.(n--l) 


'£ Cr is in Ibm/ft -hr and the pressure terms and P 
are in psi, the equation 3.6 can be written ae 



'K.n-l + J'K.n'' - l 4 t:;g- 


n-1^ 


■ G 

144 •g^ 




_.s_ 42. 

144 .g^ • T 
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